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Abstract 

A number of simplified dynamical problems is studied in an at- 
tempt to clarify some of the mechanisms leading to turbulence and 
the existing proposals to control this transition. 

A simplified set of boundary layer equations displays a solution 
that corresponds to the rolls and streaks instability and exhibits its 
streamwise localized nature. The effect of random phases as a device 
to control the transfer of energy to the small scales is studied both 
for the shell and the boundary layer models. In spite of the simplified 
nature of the models, they also provide some insight on the prospects 
for active turbulence control by external body forces. 

1 Introduction 

A realistic model for the transition to turbulence and for the mechanisms 
of turbulent energy dissipation, requires not only an accurate control of the 
fluid equations but also a consideration of the effect of different boundary 
conditions and of the external forces. Nevertheless, as in other fields of 
physics, a step towards the understanding and control of the phenomena 
is the identification of the universal mechanisms which are always present, 
independently of the particular system or the boundary conditions. One such 
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mechanism is associated to the role of coherent structures in the generation 
of turbulent energy. 

Flow visualization techniques [l]] 0, as well as two-point correlation 
measurements 0, have revealed that a large amount of the near wall turbulent 
energy is associated to coherent structures (streamwise rolls and streaks) and 
to their breakdown (bursts). Of particular importance for the shear stress 
peaks at the wall, that are associated to the bursting events, is the trans- 
fer of fast moving fluid to the near wall region (sweeps) and the transfer of 
energy to small turbulence scales at a rate faster than the normal Kolmogo- 
roff cascade. Understanding and modeling these phenomena is of practical 
importance because, stabilizing or delaying the evolution of the coherent 
structures, might decrease the energy dissipation rate. For this purpose, in 
Sect. 2, a set of simplified boundary layer equations is developed, which 
provides a good degree of analytical control over the coherent structure so- 
lutions. From the approximate equations and their solutions, the generation 
and lifting effect of the rolls-plus-streaks instability is clearly understood as 
well as its streamwise localized nature. 

Another mechanism of critical importance, in the development of tur- 
bulence, is the transfer of energy from the large to the small scales. The 
rate of transfer may be also be affected by the time scale of destabilization 
of the coherent structures. The transfer of energy between different length 
scales is favored by the establishment of precise phase relations between the 
dynamical evolution of the intervening scales. This led several authors || |J 
to suggest that randomizing the phases might have a strong effect either on 
blocking the cascade of energy transfer or on destabilizing some structures. 
In Sect. 3 the effect of random phases is studied, first in the shell model and 
then in the boundary layer equations. The study of these two cases puts 
into evidence the reasons and conditions for effectiveness of random phase 
control techniques. In the shell model, phase randomization of almost any 
mode, has an effective blocking effect on the energy transfer to the small 
scales. However, in the boundary layer equations, the effect is not so striking 
and is only observed if the phases that are randomized are those of some 
particular modes. The reason for this different behavior is the local (in mo- 
mentum space) nature of the interactions in the shell model as opposed to 
the long-range mode interactions in the boundary layer equations. The con- 
clusion is that for boundary layer turbulence control, phase randomization is 
effective only if a fine tuning of the phase-randomized modes is done. When 
the phase randomization is implemented by boundary conditions (artificial 
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roughness, for example) the need for a fine tuning of the randomized modes 
may hinder the applicability of the method. 

Finally, in Sect. 4, the prospects for active control of the coherent struc- 
tures by external body forces is considered. Here again, the fact that simpli- 
fied boundary layer equations are used, provides a clear view of the relevant 
dynamical mechanisms. Feasible control signals are derivatives at the wall 
and external forces cannot be designed with complex wall-normal dependen- 
cies. Taking these two facts into account the conclusion is that control by 
external body forces has almost no effect on the development, growth and 
break-up of the coherent structures. Therefore the turbulence production 
rate is not significantly changed. However, by controlling a small layer of 
non-active fluid near the wall, a reduction in skin friction drag might be 
obtained. 

2 Wall region coherent structures from the 
boundary layer equations 

Numerical simulations and a detailed analysis of the flow equations § have 
uncovered a possible self-sustained mechanism for the coherent structures in 
the wall region. Let x, y and z be the streamwise, the wall normal and 
the spanwise coordinates. The process starts from counter-rotating stream- 
wise rolls which would decouple from the streamwise velocity component if 
the flow were independent of x. However the rolls redistribute the momen- 
tum, creating fluctuations in the streamwise velocity (streaks). The spanwise 
inflections create an instability of the streaks and this in turn leads to a non- 
linear energizing of the original streamwise rolls. It is probable that this 
self-sustaining scheme, discussed by Waleffe|| for a sinusoidal wall-bounded 
shear flow, captures the correct physical process. However, because of its 
fragmented mathematical description, it is difficult to use it as a basis for 
quantitative predictions and, for example, for the evaluation of effective con- 
trol mechanisms. It would be better to obtain a single solution of the flow 
equations comprising all the effects (rolls, streaks and instability growth). For 
the full Navier-Stokes equations such solution does not seem easy to obtain. 
However, for the simplified boundary layer equations, the problem is not so 
hard and the mechanism may be made more explicit. In particular it will 
be seen why it is essential to take the x— dependence into account from the 



3 



start. The starting point is the Navier-Stokes equation for an incompressible 
fluid 

Or j i 

— + (US7)U = -— Vp + vAU + F (1) 

dt pm 

where F is the body force which, in the case of an ionized fluid |J, for example, 

is _ _ 

a ~ a ~ ~ 
F = —E + —U x B 

Pm CPm 



with the continuity equation 

dp r , 



Of 



+ V- (p n U)=0 (2) 



p m being the mass density, a the electric charge density and v the kinematic 
viscosity 

In orthogonal curvilinear coordinates, denote by (u, v, w) the streamwise, 
the wall-normal and the spanwise components of the physical velocity field 
U. Consider reference quantities, of the order of typical physical parameters, 
and change to non-dimensional variables 



T — -2L- 11 — JL- y — f — 7 U r ■ n, — JL- 1) — iL_ 

X Lr> y S r l Lt 2 U j' Ur ~ t 
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Typical values for the reference quantities are 

U r = 100 m s _1 , L r = 1 m, S r = 10" 3 m, p r = 1.2 Kg m" 3 , E r = 500 V 
cm~\a r = 15/iCcm- 3 , u r = 1.5xl0 -5 m 2 s -1 . Then Rl = ^ = 6.66xl0 6 

and -7^- and ^ = 10~ 6 are small quantities. 

Expressing ([I]) in the non-dimensional variables assuming the product 
kS (airfoil curvature x boundary layer width) to be small and neglecting 
terms of order R^ 1 , -fji (and ^) one obtains 

with 7 = L r^ rgr = 62.5, u; = jJ§- = 0.15 and, for the ionized fluid case, 

U r p r O r JX£ J 

F x = ^oE x ■ F y = fryaE y ; F z = ^oE z 

pm, a j-t a pm 



Split the variables into steady state (u, ...) and fluctuating components 
(«',...) 

u = u + u 

V = V + V 

/ 

w = w + w 
p = p + p 

Experimentally a sequence roll-streak-burst is a localized event in some small 
region of the boundary layer. Therefore it is reasonable to use, for the dis- 
cussion of the coherent structures in that region, a quasiparallel assumption 
for the steady-state solution, namely 

_ _ du du 

v = w = — = — = 

ox oz 

For example, for a scaling solution || 

u = u e (l - exp(-yx)) 

with X = and a = °* { 1 - r) 

By differentiation of Eqs. (3) the pressure terms may be eliminated. Then, 
keeping only the linear terms in the fluctuating fields one obtains 
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(5) 



where F i are the fluctuating components of an eventual controlling force. 
Without active control F i = 0. 

Solutions of (H) (with F' = 0) 

Let 

u = }e^ x cos{(3z)U{y,t) 



v 



ef x cos{/3z)V(y,t) (6) 



P 

with 



w = sin(f3z)W(y,t) 



U(y,t) + d -^ + W(v,t) = 
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and boundary conditions 



U(0, t) = U(5*, t) = V(0, t) = V(6*, t) 
= W(0,t) = W(6*,t)=0 

5* being a wall-normal size for the coherent structure, of the order of the 
boundary layer width. The qualitative nature of the solution does not depend 
much on the form that one assumes for the steady-state solution. In the 
following one considers an exponential profile. 

u = u e ^1 - exp (^-y-p) ) ( 7 ) 

Similar results are obtained, for example, for a linear profile 

# Zero-order solution (to be used as an initial condition for the full 
solution): 

If u is replaced by its average value u m in the boundary layer 

1 e- k 
u m = u e \ 1 - - + — 



one would obtain 



U(y,t) = c u e xt sm(^y) 

V(y,t) = -^ Cv e^(cos(fy)-l) (8) 
W(y,t) = c w e xt sm{fy) 



with c u + c v + c w = and 

(2n \ 2 

This implies that for the solution to grow in time it must decrease in x— space, 
namely 

/o \ 2 

vu I Ztt 



f (10) 

That is, the rolls-plus-streaks instability leads to a localized turbulence spot 
and x-dependence is seen to be an essential ingredient to have a sustained 
structure. The zy and zx structure of this averaged-profile solution are shown 
in Fig.la,b. 
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This approximate solution is now used as initial condition for the solu- 
tions of Eqs.(|]). Reintroducing the y-dependence of u, using the continuity 
equation and (H) and integrating once in y 

If -™0+e*(y)£ = C v (t) 

j^W - + fr(y)W = C w (t) 1 ) 

In Sect. 4, in the discussion of active control by external body forces, these 
equations will be handled in a direct way. Here however, a Fourier mode 
expansion is used, which will be convenient for the phase randomization 
discussion of Sect. 3. Let 



V(y,t) = J2 n v n (t)e 



(12) 



with e = -7f (f) 2 , X = £ and r = vu, (f )% . Then 

n d -t = (~n 3 + in^f^ + z 7 ^) v n 

+7 E ( n ' ~ ^) v n >g(n\n) + C v (t)S nfi 



n 

= {-n 2 + l k -^)w n + 1 Yw nl g{n,n) + C w {t)5 nfi 



(13) 



and 



\-k (14) 



For the t = initial condition use (i>o = — e; ^i = = |; iv\ — —w± = 
and the boundary conditions are E n = = En The evolution of 
v n and w n (for n ^ 0) is obtained from ([IB]) and -u and u> are fixed by 
the boundary conditions. The allowed time-dependence of Cy and Cw in- 
sures the consistency of the process. The absolute value of the coefficients 
in Eq. (|T2|) obtained after t = 5000 x 0.001 time steps is shown in Fig.2 for 
the parameters (7T. max = 30, 7 = 5, k = 2). Fig. 3 shows the zy— structure of 
this solution. One sees that, starting from the symmetric initial condition in 
Fig. la, as time develops, the energy of the structure moves away from the 
wall. Even Eq.@, although obtained for the averaged-profile solution, al- 
ready suggests a rate of growth increasing with the local streamwise velocity. 
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This behavior of the obtained solution matches the experimental fact that 
unstable streaks lift away from the wall. 

The solution, just studied, is sufficiently simple and well under control to 
serve as a basis for the study of control methods. In the following section 
the effect of random phases in this and other models is discussed and Sect. 4 
reports on the perspectives for active control by external forces. 

3 Random phases and energy transfer in the 
turbulence cascade 

Turbulence generation implies a transfer of energy between structures at dif- 
ferent length scales. Any mechanism that interferes with this energy transfer 
may have a controlling effect on the rate of turbulence production. A mech- 
anism of this kind was put in evidence in Ref. f|. Using a spectral code for 
a direct numerical simulation of a turbulent flow, it was found that random- 
ization of the phases of some of the Fourier modes have a strong effect on 
turbulence and leads to drag reduction. Drag reduction was also obtained 
on a recent experiment of flow over a surface with a random arrangement 
of protrusions. These effects have been interpreted as originating from the 
destruction by the phase randomization of the coherence of the large-scale 
turbulence producing structures. Interfering with the formation of rolls and 
plane waves, random phases (or random roughness) would delay turbulence 
production and bursting. 

Another important effect that is associated to phase randomization is the 
blocking of energy transfer to the small scales of turbulence. This is the 
phenomenon that will be analyzed in this section, first in the context of the 
shell model and then for the boundary layer equations. 

3.1 Random phases in the shell model 

The shell model is a toy model for the turbulence cascade. It may not be 
directly relevant for the turbulence generating phenomena in the boundary 
layer, nevertheless its study provides a good insight on the mechanisms of 
energy transfer between length scales. Let u n be the Fourier modes of a 
velocity field. Then, the G ledzer- Yamada- Ohkit ani [|I~0|j fllTf shell model is 
defined by the equation 
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J t + «n (15) 

% I k n U n+1 U n+2 - U n _ 1 U n+l - U n _ 1 U n _ 2 ) + /<5 n ,4 



The energy equation is 

, N N 
" ^ ' '2 _ ,2 I. . |2 



n=l fc=l 



To study the effect of random phases in the shell model, the numerical 
scheme used is the same as in Ref.[|l2|, namely k n = A;q2" 



■ l-e-^/3 ,. x 1 



Un (t + ^ = e -^^„(t) + — \^-g n (t) - -g(t - dt) j (17) 

with 

9n = iK (k+iK+2 - \K-iK+i - ^<-i<-aJ + f S n,4 (18) 

k = 2" 4 , v = 10- 7 , dt = 1(T 4 , / = 0.005 (1 + i) 

The system is driven at mode 4 and the other modes are energized by 
a cascade of energy transfer. The simulations start from a point in the 
attractor, the mode energy being the average over 10 4 time steps. The 
average mode energy in the attractor is shown in Fig. 4. It scales with a 
power law 

h~ 
run 



_ 2 

,2 _ h 3 



in the inertial region. 

Now, starting from a point in the attractor one lets the system evolve for 
a time At = 130 but, at each step, the phase of one or several of the modes is 
randomized. The results are shown in Figs.5a-c. For Fig. 5a different random 
phases are given to the modes 10, 11 and 12. For Fig. 5b only the phase of 
mode 11 is randomized. For Fig. 5c the same random phase is given to the 
modes 10, 11 and 12. In all cases one observes a strong blocking effect in the 
energy flow to the small scales, stronger when several phases are randomized 
but also quite large when only the phase of mode 11 is randomized. By 
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contrast if, instead of random phases, one applies a small random excitation 
to an intermediate mode almost no effect is observed. 

In this model the blocking of energy transfer by random phases has a 
simple interpretation. The energy flow through mode M is 

d N N 

- K| 2 + 2^]T^K| 2 (19) 

n=M+l k=l 

= -21m (u M u M +i (k M UM+2 + ^y^Um-i) J = 2tt m 

If um is multiplied by a random phase 

u M -> u M e i9{t) 
not correlated with Um-i, um+2 and Um+i, then 

(ttm) = 



and there is no energy flow through mode M. One also sees, from Eq.(|19l), 
that the possibility of blocking the energy transfer to the small scales depends 
critically on the local nature of the energy flow. Namely, the fact that only a 
few neighboring modes are involved in the energy transfer. That means that 
any random phase, on any mode whatsoever, will have a blocking effect. As 
it will be seen later, the situation in the boundary layer equations is quite 
different. 



3.2 Random phases in the boundary layer equations 

For the boundary layer case one evolves Eqs.(|T3D, starting from the same 
initial conditions and for the same number of time steps as in Sect. 2. Fig. 6 
shows the result when the phase of the v i component is randomized at each 
step. In this case one obtains a sizable suppression of the energy transfer to 
the higher modes. However, for the simulation parameters that were used 
(7 = 5, k — 2), randomizing the phase of mode 1 is the only case where any 
appreciable effect is observed. Randomizing the phases of the higher modes 
does not seem to block the energy transfer. 

This phenomenon may be understood from the Eqs. (|l3D . First, and unlike 
the shell model, the nature of the transfer function g(n ,n) implies that there 
are long-range effects, each mode receiving contributions from all the others. 
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On the other hand one sees from the diagonal terms in Eqs. (|T3|) that, for 
each set of parameters, there is a number of intrinsically growing modes up 
to some n max , with the other modes receiving energy from these unstable 
modes by the transfer function g(n , n). It so happens that, for the simulation 
parameters that were used, mode 1 is the unstable mode and only if one acts 
on this mode does one obtain some effect. 

In view of these two phenomena (the long-range interaction and the stable 
versus unstable mode effect) the method of random phases will be effective, 
in the control of the structures, only if it is carefully designed to act on all 
the unstable modes. It must be designed in a case-by-case basis and might 
become ineffective when the kinematical parameters change. 



4 Active control. Prospects 

Here one analyses the effect of external forces (the F' in Eq.@) on the 
growth and eventual stabilization of the coherent structures in the boundary 
layer equations. The aim is to find out whether an active control scheme 
may be developed using this type of forces. If effective control is to be 
achieved, it is the dominant dynamical effects that have to be addressed. 
Therefore it makes sense to use a setting as simple as possible to be able to 
isolate and interpret the most basic effects. In particular, one may assume 
that outside the boundary layer region there is a uniform constant pressure 
without fluctuations. Then, if F y = 0, the second equation in (f|) implies 
|^ = and therefore, from uniformity, also ff = = 0. In this case the 
pressure terms may be dropped. 

Using now the information obtained in Sect. 2, namely the fact that grow- 
ing coherent structures must be localized in the x— variable, one writes 

u'(x,y,z,t) = e~Z x -j^g(z)U(y,t) 

v'(x,y,z,t) = e-^£g(z)V(y,t) (20) 

w(x,y,z,t) = e-t x g{z)W(y,t) 

and the linearized fluctuation equations become 

W-(fty)W = F' z (21) 
-£U+%- + W = 
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Notice that the linearized equations do not fix the z— dependence. 
As before, consider an exponential profile for the steady-state flow 

u(y) = u e (1 - e- ay ) 

Eqs.(|2T|) are simple enough to even allow explicit analytical solutions. For 
example the second equation in (|2~lD reduces after a change of variable to the 
modified Bessel equation. However, better insight is obtained by displaying 
its numerical solution, starting from different sets of initial conditions. For 
example, starting from an initial condition U(y,t) ~ sin(27n/), V(y,t) ~ 
(cos(2-7n/) — 1) and W ~ sin(27ry) (Fig. 7a), as in Sect. 2, one finds once more 
that time evolution lifts this structure away from the wall (Fig. 7b). 

To control this structure would mean to damp it or at least to avoid its 
instability, growth and lifting away from the wall. Instability of the structure 
generates small scale turbulence. On the other hand when the structure lifts 
away from the wall one sees a clear growth of the V— component which is 
the mathematical equivalent of an ejection or a sweep. This exchange of 
boundary layer fluid with the external fast moving fluid increases the skin 
friction drag. 

Mathematically it would be very simple to invent forcing terms F x and 
F z which stabilize or damp the structures. For example, forces proportional 
to the second derivative would be equivalent to a change of the viscosity 
coefficient . However it is not realistic to assume that forces with the appro- 
priate y— dependence may be obtained, nor it is in general possible to detect 
the corresponding ^—variation of the fluctuating flow components. At most 
one may apply a force that is essentially constant throughout the boundary 
layer and the signal to be detected might be the value of the derivative at 
y — of the velocity fluctuations. In Fig. 7c the result of a control experi- 
ment of this type is displayed. At each time step the derivatives ^\ y =o an d 
^-\ y =o are detected and a force is applied to the fluid near the wall that is 
proportional to the difference between this derivative and the corresponding 
derivative of the initial condition. That is, the control tends to preserve the 
derivative at the wall. The range of the force that is applied is the same 
as the range of the boundary condition (a unit interval). One sees that the 
control maintains the circulation in the wall region but does not avoid the 
lifting of the structure. Actually a linear control, as used, has an overshoot- 
ing effect because, as the structure lifts away from the wall, the derivative 
naturally decreases and the control, that at first counters the sweep in the 
V— component, later causes an ejection. 
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Instead of trying to stabilize a formed roll-streak structure one might in- 
stead apply a force to kill the derivatives at the wall. In Figs.8a,b one shows 
the results of such an experiment. Starting from a small random non-zero 
initial condition in the interval y e [0, 1] , a lifted structure develops (Fig. 8a) 
which is qualitatively similar to those obtained from roll-streak initial con- 
ditions. Applying a controlling force proportional to the derivatives at y = 
, one obtains the result shown in Fig. 8b. One sees that away from the wall 
the lifted structure is similar in the controlled and non-controlled cases. The 
difference is the vanishing of the derivatives at the wall in the controlled 
case. The conclusion is that, with detection of the wall derivatives and an 
applied force in the boundary layer region, the effect over the development, 
growth and eventual break-up of the coherent structures is minimal. The 
only non-negligible effect that can be achieved is that, by suppressing the 
wall derivatives, all this activity is kept away from a small near-wall region. 
In practical terms, the conclusion is that not much effect is be expected in the 
turbulence production rate. However, by keeping a non-active small layer at 
the wall, a reduction in the skin friction drag might be obtained. The effect 
is similar to one of the interpretations of a riblet structure, namely that of 
keeping a region of non-active fluid inside the grooves. The difference is that 
here no penalty would be paid for an increased contact surface. 



5 Conclusions 

Isolating the dominating terms in the near-wall flow equations, some un- 
derstanding of the turbulence-generating structures is obtained. Also, the 
simplified equations make a clear distinction between what is important and 
what it is not, for the flow control action. 

Phase randomness (introduced by surface effects, for example) which may 
be very striking in particular models and for some kinematical regions, seems 
to lack general validity and to require fine tuning of the phases that are to 
be randomized. 

Active control by external body forces, reacting to the measurement of 
the surface derivatives, have only a limited effect on the development, growth 
and break-up of the coherent structures. Therefore the turbulence production 
rate is not significantly changed. However, by controlling a small layer of 
non-active fluid near the wall, a reduction in the skin friction drag may be 
obtained. 
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The type of externally controlled body forces that have been proposed 
in the past|§ |13| may be obtained by injecting ionized fluid (the same or 



another easier to ionize fluid) in the boundary layer and acting on it by 
electromagnetic forces. In this setting, static fields already have an effect on 
the flow profile |§. However to act on the turbulence-generating structures 
active control is required. Special attention has to be paid to the fact that, 
although the control acts on the bulk of the boundary layer, the detection of 
the signals is only realistic at the wall. Surface stress being directly related 
to wall-normal velocity derivatives, this seems to be the most appropriate 
parameter on which to base active control. Detection of field fluctuations, 
in an ionized fluid, would be an interesting possibility, but it seems a more 
remote possibility from a practical point of view. 

Figure captions 

Fig.la,b - zy (a) and zx (b) structure of the averaged-profile solution 
Fig. 2 - Fourier coefficients of the solution of Eq. ([T2|) after 5000 time steps 
(At = xO.001) 

Fig. 3 - zy— structure of the solution corresponding to the Fourier coeffi- 
cients in Fig. 2 

Fig. 4 - Average mode energy in the shell model attractor 

Fig.5a-c - Average mode energy in the shell model with random phases 
in (a) modes 10, 11 and 12, (b) mode 11 and (c) the same random phase in 
the modes 10, 11 and 12 

Fig. 6 - Fourier coefficients of the solution of Eq. (|P2"D when the phase of 
the v i component is randomized 

Fig.7a-c - Initial condition (a) and time evolved solution of Eqs.(|2~lD with- 
out (b) and with differential derivative active control (c) 

Fig.8a,b - Time evolved solution of Eqs. ([2l|) , starting from a random 
initial condition without (a) and with derivative active control (b) 
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